Probability distribution of persistent spins in a Ising chain 
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We study the probability distribution Q{n,t) of n{t), the fraction of spins unfiipped till time t, 
in a Ising chain with ferromagnetic interactions. The distribution shows a peak at n = Umax and 
in general is non-Gaussian and asymmetric in nature. However for n > Umax it shows a Gaussian 
decay. A data collapse can be obtained when Q{n,t)/L°' versus (n — nmax)L^ is plotted with 
a ~ 0.45 and 13 ^ 0.6. Interestingly, Umaxit) shows a different behaviour compared to {n{t)) = P{t), 
the persistence probability which follows the well-known behaviour P{t) ~ t^*. A quantitative 
estimate of the asymmetry and non-Gaussian nature of Q{n,i) is made by calculating its skewness 
and kurtosis. 
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In recent years a lot of work has been devoted in study- 
ing 'persistence' in dynamical systems [1]. Persistence is 
the phenomenon defined as the probability that a fluctu- 
ating non-equihbrium field has not changed its sign upto 
time t. This phenomenon has been observed in magnetic 
systems [2-6], simple diffusion [7], coarsening dynamics 
[8], various models undergoing phase separation process 
[9], fluctuating interfaces [10] etc. 

In Ising system, persistence is simply the probability 
that a spin has not changed its sign up to time t after the 
system is quenched to a low temperature from an initial 
high temperature. The fraction of the persistent spins 
Pit) here is given by 



P{t)^t-\ 



(1) 



where is a new exponent not related to any previously 
known static or dynamic exponent. In one-dimension, at 
r=0, 6 is exactly known, 61=0. 375 [2]. In higher dimen- 
sions, the persistence exponent has been obtained ap- 
proximately using both analytical and numerical meth- 
ods. 

In the numerical studies, one needs to generate dif- 
ferent random initial configurations to obtain the persis- 
tence probability P{t) which is an averaged out quan- 
tity. We deflne n{t) to be the fraction of persistent spins 
till time t which has different values for different real- 
isations of randomness such that {n{t)) — P{t), where 
() denotes average over realisations. Thus n{t) can be 
defined as a stochastic variable described by a probabil- 
ity distribution function. We have precisely studied the 
probability distribution Q{n,t) of n{t) and obtained a 
number of interesting features of the distribution for the 
one-dimensional Ising model. Probability distribution of 
random variables makes interesting studies in various sys- 
tems, e.g., for the mass of spanning cluster in percolation 
[11], random Ising and bond diluted Ashkin- Teller model 
[12], conductance of classical dilute resistor network [13], 
directed polymers and growth models [14], degree distri- 
bution in networks [15] etc. In several of these systems, 
the distribution is non-Gaussian and shows many inter- 
esting features. Certain properties like self averaging. 



multifractality etc. can be studied directly from the dis- 
tribution function [16]. Also measurements like skewness 
and kurtosis [17] from the higher moments to estimate 
quantitatively the asymmetry and departure from Gaus- 
sian behavior of the distribution are possible. 

We have considered a chain of Ising spins with nearest- 
neighbour ferromagnetic interaction and simulated it us- 
ing periodic boundary conditions. The interaction is rep- 
resented by the Hamiltonian 



H 



-^E- 



i'j+i- 



(2) 



The initial configuration is random and single spin flip 
(deterministic) Glauber dynamics has been used for sub- 
sequent updating. 
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FIG. 1. Probability distribution Q{n,t) as a function of n 
at t = 500, 1000 and 4000 for L = 700. The continuous lines 
are Gaussian fits of the form exp(— (x — xo) /a). For the fits 
shown for t = 1000 and 4000 on the right side of the peak, 
a = a,. = 1.5 X 10~^ and 1.4 x 10~^ respectively. The left side 
of the peak can be fitted to the same form with a different 
value of q: a = ai — 7.0 x 10"'' as shown in the curve for t 
=500. 
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FIG. 2. Scaled probability distribution Q{n,t)/L '^ as a 
function of (n{t) - nmax)L°-^ at t = 200 for L = 500,600,700 
and 1000 shows a collapse. 



Since a large number of configurations are required to 
get accurate data, the system sizes are restricted to 
L < 1000. Primarily, we are interested in the form of 
the distribution Q{n, t) which is plotted at different times 
against n(t) in Fig 1. We have obtained Q{n,t) for sev- 
eral values of L also. 

At all times, the distribution Q{n,t) shows a peak at 
n{t) — nmax{t)- The value of nmax{t) has negligible de- 
pendence on L. At early times, there is a Gaussian decay 
of Q{n,t) on both sides of the peak. Interestingly, the 
Gaussian behavior exp(— (a; — xo)^/q:) is followed with 
different a values on the two sides of Umax, oil on the 
left and a^ on the right. For very large i, a^ is the 
only measure possible as it is difficult to fit the function 
to a Gaussian on the left side of the peak. Usually the 
decay behaviour of distributions for rare events is of in- 
terest and we find the decay is Gaussian at all times for 
n > Umax- We observe that ctj. shows a weak dependence 
on t which becomes negligible for larger system sizes. It 



is also a function of L, ar ~ L 



ai, 



which can be 



calculated accurately for initial times, follows a similar 
scaling. In fact, the scaled distribution Q{n,t)/L" plot- 
ted against {n{t)—nmax{t))L^ with a ^ 0.45 and /3 ~ 0.6 
shows a nice data collapse (Fig 2). Even at long times, 
a fairly good data collapse can be obtained with these 
values of a and (3. 

The distribution has natural cut-offs at n = and 
n — 1. Therefore as time evolves, the distribution be- 
comes more and more asymmetric as the fraction of per- 
sistent spins decreases with time. This asymmetry is 
more apparent when the probability that there is no per- 
sistent spin, (5(0, i), begins to assume finite values. 

In fact the behaviour of (5(0, t) is quite interesting it- 
self. In Fig 3, (5(0, t) has been plotted against t/L'^ and 
the data for different L values seem to fall on the same 
curve indicating that (5(0, t) is a function of t/L^ with 
the behavior 

Q{0,t) = for t/L^ < ao 



^ for t/L"^ > ao 

where ao ~ 0.001. While t/L'^ appearing as a scaling ar- 
gument is expected, what is noticeable is the small value 
of the threshold oq. 
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FIG. 3. (5(0, t) as a function of t/L^ for four different L 
values. The solid line is a guide to the eye. 

The comparison of the behaviour of the most probable 
value Umaxit) and the average value P{t) shows consis- 
tent features. In Fig 4, we have plotted both P{t) and 
nmax{t) against t. While P{t) shows the expected power 
law decay with the known exponent 9 = 0.375, Umaxit) 
shows a different behavior, rimaxit) isl^s off faster than 
P{t) and it is not possible to fit a power law to it. 




FIG. 4. {n{t)) = P(t) (x) and nmax(t) (+) are plotted 
against time f for L = 1000. While (n(t)>(t) ~ t-O-^^s rimax{t) 
shows deviation from this behavior at finite times (the dashed 
line with slope -0.375 is fitted to (n(t))). 

One can try a power law fit only for very early times 
with a value of 6 close to 0.375. For P{t) it is known 
that the behavior P{t) ~ t~^ is valid for t < t, where 
T r^ LF' . In case the power law behaviour of rimax is valid 
for a finite time in the same sense, it appears that the 
deviation from the power law takes place at a much ear- 
lier time. We find that the so called deviation occurs at 



a value of t/L"^ ^ ao indicating that the most probable 
value fails to show persistence behaviour when (5(0, t) 
becomes non-zero. Also P{t) for i ^ oo goes to constant 
value {^ L^'^'^) while Umax goes to zero for t -^ oo. 



FIG. 5. The skewness (curve marked I), kurtosis (II) and 
scaled R„ (III) at f = 500 as a function of the system size L 
have been shown. R„/t ' (data for 4 different values of time 
have been plotted) is seen to vary as L^ ' with 9i ~ 0.51 and 
dt - 0.345. 



would imply that the distribution is flat, k shows a sat- 
uration with L for all times indicating the non-Gaussian 
behavior of the distribution. It has a positive value ^ 0.2 
to show that the distribution is peaked close to the mean. 
This value is also independent of time implying that non- 
Gaussian behavior remains constant quantitatively with 
time. In Fig 5, typical variation of s and k with L have 
been shown. 

In summary, we have obtained the distribution func- 
tion for the fraction of persistent spins n{t) for a one 
dimensional ferromagnetic Ising system. The form of the 
distribution is non-Gaussian in general. The form also 
changes with time, becoming more and more asymmet- 
ric at longer times. The most probable value rimaxit) 
shows deviation from the average {n{t)) = P{t) r^ t^^ at 
times t/L'^ > oq where ao ~ 10~^. Here we also find that 
(5(0, i) begins to take non-zero values. The system also 
shows strong self averaging. 
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One can easily calculate the higher moments from the 
probability distribution. We have studied the self aver- 
aging property and also tried to estimate the asymmetry 
and the non- Gaussian behavior by calculating the skew- 
ness and kurtosis of the distribution. 

A system is said to exhibit self averaging if Rx{L) = 
(Ax)^/(x)^ — > as i ^ oo for any physical quantity 
X [12,16] (Ax is the variance — {x^) — (x)^.) Here 
we have calculated Rn{L) — (An)^/(n(t))^ to check 
whether self-averaging is present. Our results show that 
Rn{L,t) ~ L~^' where 0; ~ 0.51 indicating strong self av- 
eraging. In fact Rn{L, t) also shows a power law increase 
with time such that i?„(i, t)/t^^ for different values of L 
shows a collapse (Fig 5). Apparently, the variance is a 
weak function of time: (An)^ ^ t^°"^. 

To measure the skewness, we calculate s{t, L) = 
M:}/{M2)^^^ where the mth centered moment is Mm — 
{{x — (x))™). s measures the asymmetry of the distri- 
bution and is zero for a symmetric distribution. Here it 
shows an increase with time as expected (e.g for L — 1000 
s ~0.18, 0.36 and 0.57 for t == 200, 500 and 2000 respec- 
tively), but shows a decay with L as well. However, the 
dependence on L weakens at longer times. We expect 
s to remain finite at i ^ oo, which is not so apparent 
from the data presumably because of the small system 
sizes considered. That there is no universal distribution 
with respect to time as in [10] is reflected by fact that s 
is time dependent. 

Kurtosis is a measure of the peakedness of the distribu- 
tion. It is studied by calculating k{t, L) = Mi/{M2Y - 3- 
For a Gaussian distribution, k{t, L) — indicating that 
the peak is at the mean value. A negative value of k 
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